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Abstract: In this review I cover recent developments concerning the construction of non-relativistic 
effective theories for perturbative heavy quark-antiquark systems and heavy quark mass def- 
initions. I then discuss next-to-next-to-leading order results on quarkonium masses and 
decay, top quark pair production near threshold and QCD sum rules for T mesons. 
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HEAVY QUARK effective field theory has pro- 
foundly influenced the way we think about 
mesons made of a heavy and a light quark Q. 
It provides a systematic expansion in powers of 
as and AQCD/mg. Perhaps more important, 
it separates hard (momenta I ^ fng) and soft 
{I ~ Aqcd) physics. Hard effects can be calcu- 
lated; spin-flavour symmetry ||^ relates soft ef- 
fects and makes the framework predictive. 

The construction of an effective theory for 
mesons made of a heavy quark and antiquark 
("onia") turned out to be more difficult. There 
are important differences between QQ systems 
and Qq systems: the former develop weak-coup- 
ling bound states in the heavy quark limit, the 
latter don't. The expansion parameter for onia 
is the velocity of the heavy quark, v, rather than 
Aqcd/wq. There is no flavour symmetry. But 
first of all, the presence of four momentum scales 
niQ, rnqv, rnqv'^ and Aqcd rather than only two 
complicates the effective theory construction. 

In a seminal paper Caswell and Lepage in- 
troduced effective field theory methods to QED 
bound state calculations The QCD equiva- 
lent, called non-relativistic QCD (NRQCD) g, 
has now become an established tool in lattice 
simulations and for describing quarkonium pro- 
duction and decay However, for perturbative 
heavy quark-antiquark systems, by which I mean 
onia that satisfy rnqv^ 3> Aqcd hi addition to 
w < 1, NRQCD is not yet optimal. NRQCD fac- 
torizes the scale mq (and needs to assume only 
that mg 3> Aqcd), but does not deal with the 



large scale hierarchy rnqv ^ rngv^. It does 
not make explicit the dominance of the static 
Coulomb force and the approximate quantum- 
mechanical nature of perturbative QQ systems. 

During the past two years progress on per- 
turbative QQ systems has been rapid. Some de- 
tails remain to be clarified, but the effective field 
theory (EFT) picture is now essentially complete. 
Several advanced applications have been worked 
out. The use of effective theory may not be com- 
pulsory for perturbative QQ systems. However, 
the gain in systematics in the expansion in v, 
transparency of language and, eventually, tech- 
nical simplification, is enormous. In this sense 
effective field theory has had as profound an im- 
pact on understanding onia as it had on under- 
standing heavy-light mesons. 

1. Effective theories 
1.1 NRQCD 

Scattering processes at momentum transfer much 
smaller than ttiq are reproduced by the non-relati- 
vistic Lagrangian H 
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+ E ^(^^rx)(x^rv^) + . . . + Aight. (1.1) 



r Q 



Loop graphs involve large momentum of order 
•niQ . The large momentum regions are accounted 
for by adapting the couplings di{h) order by or- 
der in as{'mi,). By adding more operators the 
full scattering amplitude can be reproduced to 
any accuracy in an expansion in as and v. 



The La grangian ( |l.l[) appears similar to the 
heavy quark effective Lagrangian, except that it 
contains an antiquark clone of the single heavy 
quark sector and a quark-antiquark sector (oper- 
ators made out of ip and x)? which is absent in 
HQET. Indeed, the couplings in the single quark 
(antiquark) sector are identical in HQET and 
NRQCD to all orders in Q , if the same factor- 
ization prescription is used. However, the power 
counting is generally different. For example, the 
kinetic energy term 9^/(2r7ig) is expected to be 
of the same order as id^, because non-relativistic 
QQ systems have momenta of order rnqv and 
energies of order mgf^. The situation is actu- 
ally more complicated. Because the kinetic term 
is part of the leading order Lagrangian, scatter- 
ing amplitudes computed with NRQCD depend 
on mg, V and the cut-off A in a non-trivial way. 
If one chooses A several times mqv, every Feyn- 
man diagram is a complicated function of v. This 
should be compared to HQET, where, choosing 
the cut-off several times Aqcd, every Feynman 
integral is just a number. The dependence on mg 
is fixed by the over-all power and coupling of the 
operator in the effective Lagrangian. There is 
no problem of principle with non-trivial depen- 
dence on w, if the Lagrangian is defined with a 
hard cut-off larger than raqv and if the NRQCD 
quantities are computed non-perturbatively as in 
lattice NRQCD. 

However, for perturbative calculations of per- 
turbative QQ systems this is inconvenient. First, 
one calculates too much, because one will need 
the NRQCD amplitude only to some accuracy in 
the V expansion. The full v dependence has a 
technical price. It is more difficult to compute 
an integral which is a function of v than an inte- 
gral which is just a number. Second, one would 
like to use dimensional regularization (DR). Here 
the subtleties arise, because DR gets NRQCD in- 



tegrals, written down naively, wrong. Because 
the integrand depends on mg, the dependence 
on the scale /i of DR corresponds to a cut-off 
A ^ mq instead of mg » A 3> mqv [Q. Con- 
sequently, QCD is not matched correctly onto 
NRQCD. This difficulty is related to the exis- 
tence of two scales, mqv and mqv'^ in NRQCD. 
The first attempt to separate these scales was 
made in |^ in the context of cut-off regulariza- 
tions and time-ordered perturbation theory in 
Coulomb gauge. This work introduced the im- 
portant distinction of soft and ultrasoft giuons 
(photons) and the multipole expansion for ultra- 
soft giuons, but the construction remained com- 
plicated and somewhat qualitative. Subsequent 
work ^, [ill identified different momentum 
regions that should contribute to NRQCD in- 
tegrals, but dropped the soft region identified 
in [^j. These early works introduced many of 
the important concepts that appear in a non- 
relativistic EFT construction, pointed towards 
the necessity to perform expansions of Feynman 
integrands ^, ^, |l^ in order to define NRQCD in 
DR, but did not yet provide a complete solution 
to the problem of separating the scales mqv and 
mqv^ and of defining NRQCD in DR. 

1.2 Threshold expansion 

Take an on-shell scattering amplitude of heavy 
quarks with momentum mqv and giuons with 
momentum mqv'^ in QCD. The basic problem 
is to construct term by term the expansion of 
this amplitude in v. Such a construction may in- 
troduce divergent loop integrals in intermediate 
steps, even if the original expression was finite. 
We would like to use DR to regularize these di- 
vergences and we require that any integral that 
we have to compute contributes only to a single 
order in the v expansion. Solving this problem is 
equivalent to constructing a non-relativistic ef- 
fective theory in DR with easy velocity power 
counting and all scales separated. 

The expansion method described in be- 
gins by identifying the relevant momentum re- 
gions in loop integrals, which follow from the 
singularity structure of the Feynman integrand. 
This is analogous to identifying hard, coUinear 
and soft particles in high-energy scattering of 
massless particles. For non-relativistic scatter- 
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ing of heavy quarks, one finds four momentum 
regions: 



hard (h) 
soft (s) 
potential (p) 
ultrasoft (us) 



Iq I ^ ruQ, 

Iq^I ^ mqv, 

lo ~ rnqv'^, l ^ rngv, 

lo ~ r~ rnqv^. (1.2) 



Both, heavy quarks and gluons can be hard, soft 
and potential, but only gluons can be ultrasoft. 
(In the following, "gluons" may include all other 
massless modes, i.e. light quarks and ghosts.) 

The threshold expansion is constructed by 
writing a Feynman diagram in QCD as a sum 
of terms that follow from dividing each loop mo- 
mentum integral into these four regions. The di- 
vision is done implicitly, through expansion of 
the propagators. No explicit cut-offs are needed. 
The expansion rule is that in every region one 
performs a Taylor expansion in the quantities 
which are small in that region. An immediate 
consequence of this is that every integral con- 
tributes only to a single power in the velocity 
expansion. 

To give an example of the expansion rules, 
consider the propagator of a heavy quark with 
momentum {q/2 + Iq,p + I), 



1 



ll-P - 2p- l-p'^ + qlo + y + ie 



(1.3) 



and assume that p scales as mv and y — q'^ /A ~ 
m? as mw^. When I is hard, we expand the 
terms involving p and y and the leading term 
in the expansion scales as w". When / is soft, the 
term qla is largest and the remaining ones are 
expanded. The propagator becomes static and 
scales as v~^. Notice that this means that the 
kinetic energy term in the NRQCD Lagrangian 
is treated as an interaction term in the soft re- 
gion, because it is small. When I is potential, 
the propagator takes its standard non-relativistic 
form after expansion of and scales as w^^. The 
gluon propagator takes its usual form, when the 
gluon line is soft and ultrasoft and scales as w"^ 
and v~'^ , respectively. If the gluon momentum 
is potential, one can expand Zq and the inter- 
action becomes instantaneous. If we add the 
scaling rules for the loop integration measure. 



d^l ^ (s), (p), ti*(us), we can immedi- 

ately estimate the size of the leading term from 
a given region. It is clear that these rules can 
be reformulated as an effective Lagrangian. The 
hard subgraphs give the dimensionally regular- 
ized couplings in the NRQCD Lagrangian. On 
the other hand, the distinction of soft, potential 
and ultrasoft momentum implies a manipulation 
of the NRQCD Lagrangian that is not evident 



from (1.1). 

It is instructive to discuss the physical inter- 
pretation of the terms that arise in the thresh- 
old expansion on a particular diagram. Take the 
planar, 2-loop correction to the electromagnetic 
heavy quark production vertex |p^ , i.e. the ma- 
trix element {Q{pi)Q{p2)\Ql^Q\Q)- Many of the 
possible combinations of h, s, p and us loop mo- 
mentum result in scaleless integrals, which are 
zero in dimensional regularization. Such scaleless 
integrals gain significance only in the context of 
the renormalization group, in which case one has 
to be more careful about the nature of 1/e poles. 
The non-vanishing configurations are shown in 
figure |l]. Every diagram stands for a series in v 
that arises from the expansion of the integrand in 
a particular loop momentum configuration, but 
every integral in this series contributes only to a 
single power of v. 

The h-h configuration is a 2-loop correction 
to the coefficient functions in the non-relativistic 
expansion of the current Qj^Q Ciip'^a^x + 
.... The leading term is of order a'^v^. There 
are two possibilities for having 1-loop hard sub- 
graphs. The first one (upper right) represents 
a 1-loop renormalization of the non-relativistic 
external current followed by exchange of a po- 
tential gluon. At leading order in the v expan- 
sion potential gluon exchange can be interpreted 
as an interaction through the Coulomb potential 
—asCF/r. This contribution is of order a1/v. 
The second h-p term corresponds to the insertion 
of a four-fermion operator ('0^'*^)(x^'*x) from 
the NRQCD Lagrangian. This contribution be- 
gins at order a^u. The soft subgraph in the s- 
p term can be interpreted as an instantaneous 
interaction, because a soft subgraph can be ex- 
panded in the zero components of its external 
momenta. The s-p graph corresponds to an in- 
sertion of (part of) the 1-loop corrected Coulomb 
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(h-h) (h-p) 




(h-p) (s-p) 




(P-P) (p-us) 

Figure 1: Decomposition of the planar, 2-loop ver- 
tex integral in the threshold expansion. Line coding: 
solid and curled - hard quarks and gluons, respec- 
tively; double line and wavy - soft quarks and glu- 
ons; long- and short-dashed - potential quarks and 
gluons; zigzag - ultrasoft gluons. 

potential plus higher potentials more singular in 
r. It contributes at order aj./v. The p-p term is 
dominant at small v. The double insertion of the 
Coulomb potential contributes at order a^jv'^. 
If V is counted as a^, as forced upon us by the 
dynamics of perturbative, non-relativistic bound 
states, this term is unsuppressed relative to the 
tree graph. This shows the need to partially sum 
the expansion in as to all orders. Finally, ul- 
trasoft gluon exchange is of order a^/u. This is 
actually an over-estimate. After combining all 
diagrams with an ultrasoft gluon, one finds that 
the coupling of ultrasoft gluons to heavy quarks 
has at least a factor of v. This cancellation is 
manifest in Coulomb gauge, in which only the 
spatial component of the gauge field can be ultra- 
soft. Hence the p-us term is of order o?gV. Note 
that in an expansion scheme in which u ~ Ofs, 
only four of the six terms in figure ^ are needed 
at next-to-next-to-leading (NNLO) order. 

1.3 EfTective theory again: PNRQCD 

We have now defined NRQCD and can compute 



its couplings in DR. This leads to a new set of 
graphs in which hard subgraphs have been con- 
tracted to a point. The remaining loop integrals 
can still be soft, potential or ultrasoft. One could 
introduce different fields for p, s and us quarks 
and gluons and implement the threshold expan- 
sion rules at the level of effective propagators and 
vertices But since soft heavy quarks and 

gluons and potential gluons do not appear as ex- 
ternal lines for non-relativistic systems, we can 
integrate them out. As can be seen from the 
structure of the threshold expansion, soft sub- 
graphs have the same conibinatorical structure as 
hard subgraphs. They can be contracted to effec- 
tive operators, which are non-local in space, but 
local in time (instantaneous) . These effective op- 
erators provide a generalized notion of the heavy 
quark potential. The resulting effective theory 
contains only potential quarks and ultrasoft glu- 
ons. Velocity power counting is then trivial. The 
scheme is as follows: 

-Cqcd \Qih, s,p), g{h, s,p, us)] fi > m 

i 

-Cnrqcd [Q{s,p), g{s,p, us)] mv < iJ < m 
I 

-CpNRQCD [Q{p), 9{us)] fi<mv 

Such a construction has been proposed first, by 
tree level matching, in I follow in calling 
the second EFT potential NRQCD (PNRQCD). 
In the context of the threshold expansion, which 
provides a matching prescription for loop graphs, 
PNRQCD has been discussed in [ |l5[ [l6| . A some- 
what different, but probably conceptually equiv- 
alent construction has been proposed recently in 

0- 

In general the PNRQCD Lagrangian can be 
written as 

•^PNRQCD — -^NRQCD + -^non-local- (1-4) 

'Cnon-iocal collccts all non-local interactions. The 
local interactions are exactly those of NRQCD, 
but the interpretation is different, because only 
potential heavy quarks and ultrasoft gluons are 
left over. In diagrams constructed from >Cnrqcd 
the gluon propagators are always expanded ac- 
cording to their ultrasoft scaling rule and the 
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heavy quark propagator has the famihar non- 
relativistic form, while in diagrams constructed 
from £nrqcd gluons are also soft and poten- 
tial and the heavy quark propagator can also 
be static. Because only potential quarks and ul- 
trasoft gluons are left in PNRQCD, the interac- 
tion terms have definite velocity scaling rules. A 
potential quark propagator in coordinate space 
scales as v^, so a quark field in PNRQCD scales 
as u'^/^. An ultrasoft gluon field counts v'^. 

Integrating out potential gluon exchange be- 
tween a quark and an antiquark at tree level gives 
the leading order Coulomb potential. The unper- 
turbed PNRQCD Lagrangian is 



•'-PNRQCD 



(1.5) 



2mr 



2mc 



Since v ~ as^rngv) is assumed, all terms scale as 
v^; the Coulomb interaction cannot be treated as 
a perturbation as is of course anticipated. One 
can rewrite the PNRQCD Lagrangian in terms 
of a 'tensor field' [ip (g) x^{t^R^''^) that depends 
on the cms and relative coordinates. The un- 
perturbed Lagrangian describes free propagation 
(with mass 2m) in the cms coordinate. The prop- 
agation of [ip (g) x^](^j ^) ill its relative coordi- 
nate is given by the Coulomb Green function of 
a particle with reduced mass m/2. In calculat- 
ing diagrams with Coulomb Green functions, one 
sums corrections of order (ag/i))" to all orders. 
The remaining terms can be treated as pertur- 
bations in V and as around the unperturbed La- 
grangian. 

Higher order non-local interactions follow af- 
ter matching potential gluon exchange to bet- 
ter accuracy and after integrating out soft loops. 
We can match assuming ^ u, and treat the 
Coulomb potential as a perturbation when match- 
ing NRQCD on PNRQCD. A general non-local 
operator is a function of r and derivatives di act- 
ing on the (anti)quark field. It is non-polynomial 
in r, but polynomial in di. In general, it may 
have ultrasoft gluon fields attached to it. Non- 
local operators are singular, for example (r) 
or aj./r'^. These singularities are harmless, be- 




tWVWV 



(a) 



SV{r) 




(b) 



Figure 2: (a) s-s-s region that gives rise to an in- 
frared logarithm (left); p-p-us region which contains 
the corresponding ultraviolet logarithm, (b) In PN- 
RQCD notation the two NRQCD graphs in (a) are 
interpreted as a (Coulomb) potential insertion (left) 
and an ultrasoft 1-loop diagram (right). The shaded 
bar represents the propagation of the QQ according 
to the Coulomb Green function. Line coding as in 
figure |l|. 

cause PNRQCD is defined with a cut-off (we im- 
ply dimensional regularization) and more singu- 
lar non-local operators are in fact suppressed in 
V. The dimcnsionally regulated quark- antiquark 
potential up to order w^, projected on a colour 
singlet, spin-1 QQ pair can be found in p^ . 
Potentials are short-distance coefficients, i.e. 



they are for PNRQCD what the di in (1.1) are 
for NRQCD. Consider the Coulomb potential in 
(1.5). Its coefficient is renormalized such that 



as — > as{n)vc{as, ^J.r). 



(1.6) 



The coefficient function Vc{as, ^ir) contains loga- 
rithms of /xr, which are analogous to logarithms 
of the heavy quark mass in the coefficient func- 
tions of the NRQCD Lagrangian. Up to order 
a^, to which Vc{asT^Jir) is known exactly [ p8| , all 
logarithms In /ir can be absorbed into as{\/r). 
This is no longer true at three loops. There ex- 
ists an uncancelled infrared divergence in the left 
diagram of figure |^ ||l^ which, after subtrac- 
tion of the pole in DR, gives rise to a logarithm 
not related to the running coupling. This loga- 
rithm is analogous to a non- vanishing anomalous 
dimension of local operators in the NRQCD La- 
grangian. The scale dependence is cancelled for 
a physical process by the scale dependent PN- 
RQCD matrix element, in this particular case 
the right diagram of figure §1. In PNRQCD 



5 



Heavy Flavours 8, Southampton, UK, 1999 



M. Beneke 



VWvWv 



Figure 3: NRQCD graph that generates the mixed 



non-local/uhrasoft interaction in (1.8). Line coding 
as in figure |l| 

notation, which docs not make use of potential 
and soft lines, the correspondence is shown in 
figure ||b. The interpretation of potentials as 
matching coefficients implies that the Coulomb 
potential is not identical to the static potential 
defined as the vacuum expectation of a Wilson 
loop in the limit T ^ oo . The Coulomb po- 
tential is logarithmically sensitive to the ultrasoft 
scale, but it is not infrared divergent. This state- 
ment applies to any other potential. 



Although the local interactions in (1.4) are 
the same as those in NRQCD, we must write 
them in a different form to account for the ex- 
pansion rules for ultrasoft gluons. When an ul- 
trasoft gluon line with momentum I connects to 
a quark line with momentum k — 1/2 for the in- 
coming and k + 1/2 for the outgoing quark line, 
the threshold expansion instructs us to expand 
the quark-gluon vertex and quark propagator in 
l/k V. All gluon interaction terms in 'CJi^pQ^Q 
should be understood as multipole-expanded, for 
instance |||, |l^ 

[tP^A'd'iP] {x) = ij\x)A\t,Q)d'ij{x) 

+ i:\x) {x^d^A'{t, 0))d''4!{x) + . . . , (1.7) 

and likewise for all other terms in the NRQCD 
Lagrangian. The leading ultrasoft interactions 
follow from multipole expansion of the gauge field 



terms in ip"^ {iD'^ + D'^ /{2mQ))tfj in (1.1) together 
with a non-abelian non-local term that comes 
from matching the graph in figure|3|on PNRQCD. 
The following collects all ultrasoft interactions up 
to order v^^^'^: 

Cus - 9s[4>^T^^]{x)A''^{t,Q) 

+ gs['i}^^T^M{x)x'd'A''^{t, 0) 



- 9')r^^](a;)A*^(i,0) 



2mQ 
antiquark terms 



•(-^)5./''^^r-M^^(i,0) (1.8) 

The first line scales as v^^^ relative to the lead- 
ing order terms in the PNRQCD Lagrangian. 
The other three interaction terms scale as u'^/^. 
Using [x^,d'^] — —29* and the equation of mo- 
tion at leading order in u, which includes the 
Coulomb potential, the w^^^ interactions com- 
bine into a chromo-electric dipole operator up to 
higher order terms. (Note that this shows that 
the distinction of non-local and local operators 
in PNRQCD is ambiguous, because they can be 
converted into each other by the equation of mo- 
tion. Likewise a classification in powers of 1/toq 
is not useful.) Introducing the ultrasoft covariant 
derivative D'^ = 9" — igsA'^{t, 0), we obtain 



^PNRQCD = 

4>\x) iiD" + 



(1.9) 



2mr 



gsx'E\t,Q) i^{x) 



antiquark terms 



d?f [VlVfc] (x + r) VabMr. d') [xlxd] {x) 



where 



(1.10) 

This Lagrangian is manifestly invariant under ul- 
trasoft gauge transformations U{t,Q). Note that 
the spatial components of the gauge field trans- 
form covariantly under ultrasoft gauge transfor- 
mations. Beyond tree-level the coefficient of the 
chromo-electric dipole operator receives correc- 
tions that can be computed in an expansion in 
Us- The PNRQCD Lagrangian ( |1.9| ) appears to 
be equivalent to the Lagrangian derived in [l4| . 
The difference is that introduces a colour de- 
composition of the tensor field [ip ® x^](^:-^7^) 
and expresses the Lagrangian in terms of a colour 
singlet and a colour octet field 5" and O. 

A Green function with no external ultrasoft 
lines requires at least two ultrasoft interactions. 
The leading ultrasoft correction is order v from 
the ultrasoft covariant derivative in (^). How- 
ever, we can use an ultrasoft gauge transforma- 
tion to gauge A" away, and hence the leading 
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ultrasoft correction is order v^. The Lagrangian 



(|l.9| ) can be used to compute all leading ultrasoft 
contributions. In particular, ultraviolet renor- 
malization of ultrasoft graphs cancels the scale 
dependence of the potentials to order v^. In |Q 
the scale dependence of the potentials has been 
computed using this correspondence. 



1.4 Non-relativistic renormalization group 

The scale hierarchy mq 3> rnqv 3> tiiqv^ also 
implies large logarithmic corrections Us In v and 
aglnv^ to the coefficient functions. These loga- 
rithms can be summed to all orders in perturba- 
tion theory as follows: 

1) Match QCD and NRQCD at a scale fih ~ 
TOQ, i.e. compute the coefficients di as ex- 
pansions in as{fih)- 

2) Compute the renormalization group scaling 
of the di in NRQCD and evolve them to a 
scale /is ~ rnqv ^ 1/r. 

3) Match NRQCD and PNRQCD at the scale 
fj,s, i.e. compute the potential coefficients 
Vi as expansions in ^^(Ms) and in terms of 
the di{^s)- 

4) Compute the renormalization group scaling 
of the potentials in PNRQCD and evolve 
them to a scale ~ mgv^ ~ rnQO^. 

5) Compute the PNRQCD matrix elements 
with ultraviolet subtraction scale ^us- 

I briefly discuss items 2) and 4), but note that 
an explicit calculation remains yet to be done. 

NRQCD contains a single-heavy quark sec- 
tor, which is identical to heavy quark effective 
theory (HQET). Since heavy quark-antiquark op- 
erators do not mix into this sector, its renormal- 
ization can be discussed separately. Operator 
renormalization in NRQCD arises from ultravio- 
let divergences in potential and soft loops. The 
single-heavy quark sector has no potential loops 
(all quark poles on one side of the real axis in 
the complex plane of loop momentum zero com- 
ponents); the anomalous dimension matrices are 
identical to those in HQET. It is convenient not 
to introduce non-local time-ordered product op- 
erators as is usually done in HQET, but to have 



lower dimensional operators mix into higher di- 
mensional ones. For example, a vertex diagram 
with two insertions of the chromo- magnetic oper- 



ator of (1.1) requires a counterterm proportional 



to the Darwin interaction, hence 



Mil — d2 



bag 
"2^ 



di 



(1.11) 



etc.. The single- heavy quark sector mixes into 
heavy quark-antiquark operators. These opera- 
tors are renormalized by soft and potential loops. 
Mixing through potential loops is responsible for 
the scale dependence of the quark-antiquark cur- 
rent ip'^a^x that appears first at two loops ||2^, 
p3| . Higher dimension operators in the single- 
heavy quark sector can mix into lower dimension 
operators in the heavy quark-antiquark sector 
through potential loops, because potential gluon 
exchange can contribute factors mQ/\p\, where p 
is a relative momentum of order rnqv. This never 
happens for soft loops. One can introduce two 
separate renormalization scales fix and ^2 for soft 
and potential loops and compute the correspond- 
ing anomalous dimension matrices. NRQCD co- 
efficient functions Ci{^i, ^2) then depend on both 
scales. In the end we only need Ci{^,fj.), which 
evolves with the sum of the two anomalous di- 
mension matrices and the distinction is not nec- 
essary. However, the fact that potential loops 
are infrared finite tells us that the evolution in 
/j,2 stops at a scale of order mqv. Soft loops are 
not always infrared finite; this connects the evo- 
lution in fii to the ultrasoft region. 

Identifying = ^2 the NRQCD evolution 
terminates at the scale /is at which one matches 
to PNRQCD. The evolution of the potentials and 
other operators such as the chromo-electric dipole 
operator is then determined by the ultraviolet be- 
haviour of ultrasoft loops. 

A different implementation of renormaliza- 
tion group scaling from the one presented here 
has been suggested in ||l^. In this case, too, 
an explicit calculation of operator renormaliza- 
tion has not yet been performed and the equiva- 
lence of the two approaches remains to be demon- 
strated. 

1.5 Aqcd ^ rnqv'^ or larger 

The theoretical framework described so far re- 
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quires mQV^ ^ Aqcd- Nice as it may be, it 
can be applied safely only to extremely elusive 
systems such as toponium. There will be non- 
perturbative corrections suppressed by powers of 
Aqcd/("t.qw^) in addition to what I discussed 
above, which can be estimated by applying an 
operator product expansion to PNRQCD matrix 
elements ^J]. In case of toponium these correc- 
tions are estimated to be small for the inclusive 
production cross section. Note that there are no 
non-perturbative modifications of the potentials, 
because they are short-distance objects. 

What happens if mgw^ is not that large? 

If niQ ^ Aqcd ^ rngv, PNRQCD makes no 
sense. One can use NRQCD, but the NRQCD 
matrix elements are non-perturbative. 

If niQV ^ Aqcd ^ TTT-Qv'^, one can match 
perturbatively (in as{mQv)) to PNRQCD, but 
since as{mv'^) ~ 1 the self-coupling of ultrasoft 
gluons is unsuppressed. Velocity power count- 
ing is different from the one above which used 
as{mv'^) ~ as{mQv) ^ v. Non-perturbative glu- 
ons screen the ultrasoft scale mqv^ and Aqcd 
takes its role. The coupling of ultrasoft gluons 
to heavy quarks is non-perturbative but small, of 
order Aqcd/(™qw)- Hence ultrasoft effects en- 
ter a heavy quark-antiquark scattering amplitude 
as an uncalculable non-perturbative contribution 
beginning at order (Aqcd/(™qw))^. Up to this 
accuracy the amplitude is still determined by po- 
tentials. In particular, the potentials and the 
scale dependence of PNRQCD matrix elements 
remain perturbatively calculable. 

2. Heavy quark mass definitions 

Before turning to applications of non-relativistic 
field theory, I discuss the concept of the heavy 
quark mass. The (P)NRQCD Lagrangian is nor- 
mally expressed in terms of the quark pole mass 
and this has been assumed so far. If not I should 
have added a term SmQip''ip to the (P)NRQCD 
Lagrangian. There are good reasons to make use 
of this option, related to the fact that the pole 
mass, though an infrared safe quantity in per- 
turbation theory [^, incorporates uncalculable 
long-distance contributions of order Aqcd- Orig- 
inally discovered and discussed in the context of 
HQET m, 1^, Hi, the problem is acute when- 



ever heavy quarks are not off-shell by an amount 
of order TOq. Recently there has been renewed 
interest in this problem in the context of QQ sys- 
tems |2^, |3^ and suitable alternative definitions 
of the heavy mass have been proposed and used 
in applications. Such masses have the generic 
property that they differ from the pole mass by 
an amount linear in a subtraction scale fij |28). 
In this section I review the currently used mass 
definitions. 

2.1 MS and pole mass 

The MS mass mg (/i) is the coefficient of the oper- 
ator QQ in QCD subtracted in the MS scheme. It 
is best understood as a coupling constant just as 
as- The MS mass is scale-dependent. The scale 
dependence is related to loop momenta I ^ mg. 
For this reason, it makes no sense to evolve the 
MS mass to scales parametrically smaller than 
niQ. While formally possible, this generates fake 
logarithms of the ratio of scales. However, mg = 
rnQijriQ) is a very useful reference parameter, 
just as as{mz)- 

The pole mass is the location of the pole in 
the full heavy quark propagator. (The weak in- 
teractions are switched off, so that quarks are 
stable. A finite decay width would not alter the 
conclusion of this subsection [0.) As such it is 
defined order by order in perturbation theory. Its 
relation to the MS mass is given by 

^ = 1 + Vfc„a,(mQr. (2.1) 

For h quarks, neglecting internal the charm quark 
mass effects, 

fci = 0.424, fca = 0.940, 

fcg = 3.096, fc4 = 13.60 (2.2) 

The first two coefficients are known analytically 
HI], the third is known numerically |Q. For 
the fourth order coefficient I have used the so- 
called "large-/3o" estimate |Q , which turned out 
to approximate fc2,3 very well. 

Neither of the two masses is a useful pa- 
rameter for perturbative calculations of QQ sys- 
tems. The MS mass would imply Smq ^ mqag 
and hence Smqip'^ip ^ would dominate the 
(P)NRQCD Lagrangian. One must have Smq ~ 
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niQV^ or smaller. The pole mass is long-distance 
sensitive at order Aqcd • There would be nothing 
wrong with this if not for the following two facts: 
first, the quantities we would like to compute 
with non-relativistic field theory such as correla- 
tion functions of two heavy quark currents near 
threshold are less long-distance sensitive than the 
pole mass. This desirable property is lost, if one 
uses renormalization conventions which are more 
long-distance sensitive than the quantity of inter- 
est. Second, long-distance sensitivity is related to 
large perturbative corrections through infrared 
renormalons js^. Here this means that series 
expansions diverge as {2f]o)"'nln^ , where /Sq is de- 
fined through das/d In ^"^ — —Poo^l + . . .. The di- 
vergence enters less long-distance sensitive quan- 
tities only through fc„+i ~ {2Po)"-nln^ and would 
be much milder, if a suitable mass renormaliza- 
tion convention were implied. For the following 
discussion it is useful to introduce an "asymp- 
totic counting" for perturbative coefhcients. The 
asymptotic counting of the coefficients (^^) is 
kn ~ (n — 1)! fi/niQ, where the definition of kn 
has been generalized to an expansion in as{i-i) 
and we neglect the factors (2/3o)" and in this 
schematic notation. I will also call a series ex- 
pansion "convergent" , if it diverges less rapidly 
than the kn- Coefficients of "convergent" series 
count as order 1 and terms in a convergent series 
are counted only according to their power in as- 
Asymptotic counting may appear abstract 
and irrelevant to next-to-next-to-leading order 
calculations. The physical systems which we dis- 
cuss later (and numerous quantities related to B 
meson decays) show that this is not so. This has 
led to heavy quark mass definitions, which satisfy 
SniQ ~ mQv'^ and are convergent in asymptotic 
counting. 

2.2 PS mass 

The potential subtraction (PS) scheme |Q is 
based on the observation that there is a cancel- 
lation of divergent series behaviour in the combi- 
nation 2mQ + [l^(r)] Coulomb- This can be seen ex- 
plicitly at 1-loop and in the large-/3o approxima- 
tion ^ by combining the results of ^ 
and by a diagrammatic argument at two loops 
||29|| and, probably, in higher orders. The poten- 
tial subtracted (PS) mass at subtraction scale fif 



is defined by 

TOQ,Ps(M/) = "IQ + i J -^-^ [V'(g)]couiomb 
= niQ - fif^ln{fif/fi)as{iJ-y'^^, (2.3) 



n=0 



where 



[^^(9)]couiomb = ^Vc{as,ql n) (2.4) 



is the Coulomb potential in momentum space de- 
fined as a PNRQCD coefficient function as dis- 
cussed above. The coefficients /o,i,2 are given in 
|2£]. The large-/3o estimate is easily derived from 
1 3£ ] and I will use the result for l-^ below. Because 
the Coulomb potential is scale-dependent at or- 
der and beyond, it depends on the PNRQCD 
matching scale. The coefficients /„, n > 2, inherit 
this scale, which has to be specified in addition 
to fif. We can choose this scale equal to /i/, so 
that the logarithm in the Coulomb potential |2^ 
does not contribute to the PS mass. 

In asymptotic counting Z„ ^ n\ ^i/nf, so that 
the combination rnqkn — fifln-i, which appears 
in the relation of mg^ps to tTiq, is convergent as 
desired. This is true for any ^J,f > few x Aqcd- 
To satisfy Smq mqv^ or smaller, must not 
be parametrically larger than raqv. It is useful to 
choose fJ.f oi order mqv ~ niqas- Note that this 
implies that mqkn and fifln~i are of different 
order in as, but cancel asymptotically. 

2.3 IS mass 

The IS scheme was proposed in and uses 
one half of the perturbative T(15) mass as quark 
mass parameter. The perturbative T(15') mass 
is related to the physical T(15) mass Afx(is) by 



M- 



T(1S) 



2mr 



.IS 



A. 



(2.5) 



where Ay is a poorly known non-perturbative 
contribution, which is most likely less (if not con- 
siderably less) than 100-150 MeV. (Some versions 
of the IS scheme eliminate the bottom quark 
mass in favour of the physical T(IS') mass, the 
advantage being that the input parameter is a 
physical quantity which is very accurately mea- 
sured. I prefer the version stated above, because 
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it does not obscure the presence of an unknown 
non-perturbative contribution. The IS scheme 
can also be defined for top quarks |3^. It uses 
the perturbative toponium IS mass under the as- 
sumption of a stable top.) Parametrically Ax is 
of order {Aqcd / (nT-Qv))^ and not of order Aqcd- 
This guarantees that mg ig is less long-distance 
sensitive than the pole mass for perturbative QQ 
systems. 

The IS scheme does not have an explicit sub- 
traction scale fi f , but it is similar (up to "finite" 
renormalizations) to the PS scheme with /i/ of 
order niQas- The series expansion is 



mQ.is/mQ = 1 - ^ e„(^) Q!s(m) 



n-l-l 



(2.6) 



n=0 



The coefficients are known exactly up to 62 |Q 
and will be given in Sect. below. The asymp- 
totic counting is e„ ~ n\ fi/ (rngas), where the 
factor iJi/{mQas) follows from the fact that the 
physical scale is of order mqas- Hence the com- 
bination fc„ — e,i_iQ!s, which enters the relation of 
'7^Q,ls and mg as coefficient at order a" forms 
a convergent series with coefficients of order 1. 
As expected from the correspondence to /i/ ~ 
mqas, one has to combine coefficients of differ- 
ent order in as [ p7[ . 

Strictly speaking the IS scheme cannot be 
consistently used in NNLO (that is, keeping 620;^ 
in ( |2.6| )) and beyond for bottom quarks, since the 
T(IS') ultrasoft scale mgal is of order Aqcd- As 
discussed in Sect. in this case the leading ul- 
trasoft contribution (which would be order for 
a perturbative system) is non-perturbative and 
of the same parametric order as the term 620;^. 
This should be kept in mind, but in the following 



I use (2^) as a formal definition of the scheme, 
including the NNLO term 62- 

2.4 Kinetic mass 

The so-called kinetic mass has its roots in B 
physics The B meson mass has the heavy 

quark expansion 



TTiB = rrib + A 



(2.7) 



with /i^ and fiQ related to the matrix elements 
of the kinetic energy and chromo- magnetic oper- 
ators, respectively. The kinetic mass can be un- 



derstood as a perturbative evaluation of this for- 
mula, in which the matrix elements include loop 
momentum integration regions below the scale 

TOQ,kin(M/) 

>^(M/)" 



rriQ 



[A(M/)] 



pert 



2m, 



Q 



(2.8) 



pert 



The matrix elements on the right hand side sub- 
tract the long-distance sensitive contributions to 
the pole mass order by order in iif/mq and a^. 
While easily stated, the definition of power di- 
vergent matrix elements in perturbation theory 
is largely arbitrary. The convention for the ki- 
netic mass used in the literature uses an indirect 
definition through heavy flavour sum rules ]40| , 
which is rather complicated when compared with 
the other two mass definitions above. The rela- 
tion between the kinetic and the pole mass is 
known exactly at order (including terms of 
order ^"j/niQ) and in the large- /3o limit at order 

2.5 Comparison 

In table |l] I compare the various mass defini- 
tions for b quarks numerically using the MS mass 
as a reference parameter. Each entry gives the 
value of the mass using a l-loop/2-loop/3-loop/4- 
loop relation. For reasons that will become clear 
later it is interesting to have four-loop accuracy. 
Where available I have used large-/3o estimates 
to obtain the four-loop value. I used as{mt) as 
perturbative expansion parameter (with one ex- 
ception, see below). In defining the kinetic mass, 
terms of order {fif/nib)^ or less are dropped. At 
order the large-/?o estimate is used, dropping 
all known terms that are subleading in this limit. 
(The relevant formula is given in the preprint ver- 
sion of @.) 

The IS mass is computed in two different 
ways. First, I express it as a series in as{rnh) 
with coefficients fc„ — en-ictsirnb) as explained 
above. A 4- loop relation then requires 63, which 
is not yet known, or at least its large-/3o value. 
This result is shown in brackets in table |l|. The 
second way first computes the PS mass at the 
scale fif = 2GeV as series in as{Wib) and then 
relates the IS mass to the PS mass as an expan- 
sion in as(2 GeV). In this case 63 is not needed. 
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rrib 


mb,ps(2GeV) 


m-h.kin(l GeV) 


mbss 


^bjpolc 


as{mz) = 0.118 [q44.25 GeV) = 0.2240] 


4.15 


4.36/4.44/4.47/4.48 


4.41/4.49/4.50/- 


4.36(50) /4.60(62)/4.67(66)/4.66(-) 


4.55/4.75/4.89/5.04 


4.20 


4.41/4.49/4.52/4.54 


4.46/4.54/4.56/- 


4.41(55)/4.66(68)/4.72(72)/4.7i(-) 


4.60/4.80/4.95/5.09 


4.25 


4.46/4.55/4.58/4.59 


4.52/4.60/4.6i/- 


4.46(61)/4.71(73)/4.78(77)/4. 76(-) 


4.65/4.85/5.00/5.i5 


4.30 


4.52/4.60/4.64/4.65 


4.57/4.65/4.67/- 


4.52(66) /4.76(78)/4.83(82)/4.82(-) 


4.71/4.91/5.06/5.20 


4.35 


4.57/4.66/4.69/4.70 


4.62/4.71/4.72/- 


4.57(71)/4.82(84)/4.88(88)/4.87(-) 


4.76/4.96/5.11/5.25 


as{mz) = 0.121 [as(4.25GeV) = 0.2355] 


4.15 


4.37/4.45/4.49/4.5i 


4.42/4.51/4.52/- 


4.37(52)/4.63(65)/4.70(69)/4.68(-) 


4.57/4.79/4.96/5.i4 


4.20 


4.42/4.51/4.55/4.56 


4.48/4.56/4.58/- 


4.42(57) /4.68(70)/4.75(75)/4.73(-) 


4.62/4.84/5.01/5.i9 


4.25 


4.47/4.56/4.61/4.62 


4.53/4.62/4.64/- 


4.47(62)/4.73(76)/4.81(80)/4. 79(-) 


4.67/4.90/5.07/5.25 


4.30 


4.53/4.62/4.66/4.67 


4.58/4.67/4.69/- 


4.53(68) /4.79(81)/4.86(85)/4.84(-) 


4.73/4.95/5.12/5.30 


4.35 


4.58/4.68/4.72/4.73 


4.64/4.73/4.75/- 


4.58(73) /4.84(86)/4.92(91)/4.89(-) 


4.78/5.00/5.18/5.35 



Table 1: Comparison of b quark masses for given MS mass mj, at the scale mt, for two values of asimz)- We 
used Uf = 4 and put rric — 0. Slanted numbers make use of large-/3o estimates. All numbers in GeV. 



The result is shown without brackets in table |l|. 
The two ways of computing mb,is correspond to 
the mass analyses performed in |42j and |Q , re- 
spectively. 

Except for the b quark pole mass all other 
masses given in the table are related to the MS 
mass by "convergent" series. The difference is 
clearly visible in the magnitude of successive per- 
turbative terms. For the pole mass the minimal 
term is reached around 3-4 loops suggesting that 
its best accuracy cannot be smaller than about 
150 MeV. On the contrary, the perturbative ex- 
pansions of the other masses behave extremely 
well. At fourth order the error in relating them 
to mfc is of order 10 MeV! 

At this level, a comment on the effect of keep- 
ing the mass of internal charm quarks is nec- 
essary. Charm quark mass effects at order 
are known for the pole quark mass [Q. For 
given rrib = 4.25 GeV and rUc = (1.1 - 1.4) GeV 
the pole mass increases by about 8-10 MeV. The 
charm mass effect on the PS and IS mass is easily 
computed in terms of the charm quark contribu- 
tion to the photon vacuum polarization. I find, 
again for given nib, that m;,_ps(2 GeV) and rrib^is 
are reduced by 11-12 MeV and 4-5 MeV respec- 
tively. These corrections could be applied to ta- 
ble I], which is generated for nic = 0. Charm 
effects at order and beyond are not known 
and introduce an error of perhaps 10 MeV. If the 
average internal loop momentum becomes small, 



the charm quark decouples and one can switch 
to an effective description with only three light 
flavours Q. 

2.6 Application: inclusive heavy quark de- 
cay 



Before returning to QQ systems it is interesting 
to see how these mass definitions fare in decays of 
B mesons. As an example consider the inclusive 
semi-leptonic B decay B Xulv. Many more 
examples can be found in | |3^ , though restricted 
to the IS scheme. (An earlier 1-loop analysis 
with the kinetic mass was done in p4[.) 

It is known that inclusive B decays are less 
long-distance sensitive than the pole mass ||2^ , 
^ and so another mass parameter is warranted. 
But the situation is different from onium sys- 
tems, because there is no scale other than urib. 
Hence mb{rrTb) is a legitimate choice. However, 
this works well only asymptotically but fails in 
low orders in as The decay rate, neglecting 
power corrections, is 



Gj,\Vub\^M^ 



1927r3 



1-61-62 + 0{al) . (2.9) 



With the second order correction in the pole mass 
scheme from the 1-loop and 2-loop correc- 
tion in various other mass renormalization con- 
ventions is computed and shown in table |[ All 
"alternative" mass definitions introduced above 
do very well. 
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M 


Si 


S2 




+0.17 


+0.10 


nib 


-0.30 


-0.13 


m6,Ps(2GeV) 


-0.03 


-0.01 




+0.11 


+0.03 


m6,kin(l GeV) 


+0.03 


-0.005 



Table 2: First and second order perturbative terms 
in inclusive b ^ u decay with various mass parame- 
ters. 



It is perhaps not evident why this is the case 
given that there is no a priori reason not to use 
the MS mass. A plausible argument is that while 
the characteristic scale in inclusive 6 — )■ u decay 
is parametrically of order TO6, it is numerically 
smaller [Q. This forces the b quark closer to 
its mass shell although not close enough to jus- 
tify the use of the pole mass. Similar improve- 
ments by using "alternative" mass renormaliza- 
tion schemes occur for other heavy quark decays 
l37i. 



3. Quarkonium 

In this section I review results on quarkonium 
bound states. Unfortunately, none of the ob- 
served charmonium and bottomonium states is 
truly perturbative, i.e. satisfies mqv^ 3> Aqcd- 
The best candidate to try our theory are the 
T(15) and r/b states. Higher excitations are al- 
most certainly non-pcrturbative, although more 
non-relativistic. 

3.1 Masses 

The quarkonium binding energy is of order mqa^ 
in leading order. For an arbitrary QQ[nl] state 
the energy is known to order For nS 

states the result of |3^ has been confirmed by 
H, |0[ |l|l. The T(15) mass, expressed in terms 
of the b quark pole mass, is given by 



4: 9 9 



1 + — 

IT 



625 



203 ' 

'Uy ' 7r2 V 48 
24537r2 1235C(3) 



216 



36 



1429 
24 

7211 

' To^ 



25 
'32 



2mb - - mlal 



1 + 1.08 + 1.76 



(3.1) 



where I = ln(16m^Q!^/(9^^)) and ag = as(^) 
and the second line is given for ^ such that / = 
(for wLb = 5 GeV), in which case as(^) = 0.30. 
Corrections to this result are order and order 



Aqcd 



Aqcd A 
mbol ) 



2+2n 



(3.2) 



from the operator product expansion 



An interesting application of (3.1) is to use 
it to determine the b quark mass. The non- 



convergence of the series ( 5T ) seems to make this 
impossible. However, if the series is expressed in 
terms of the PS mass or the MS mass, the con- 
vergence is dramatically improved as seen from 
the column referring to mb,is in table |^. In |^ 
this has been used to obtain 

mbirfib) = (4.24 ± 0.09) GeV. (3.3) 



The central value is obtained by varying /i from 
1.25 GeV to 4 GeV and by symmetrizing the er- 
ror. The total error is dominated by the unknown 
non-perturbative contribution from ultrasoft glu- 
ons and the fact that the OPE series in n of (3.2) 
may not converge. Note that (3.3) uses the 4- 



loop relation in table |l] because the series (|3.lD 
in terms of the PS or MS mass is convergent in 
asymptotic counting and hence determines the 
quark mass to order a^. This is the main reason 



why (3.3) differs substantially from the value ob- 
tained in , where nib is determined via the b 
pole mass and by a 2-loop relation. 

There exist partial results at order which 
may be used to estimate the perturbative error 
on Mx(is). The mass correction from ultrasoft 
gluon exchange has been obtained in ||51[: 



SM^\^g^^ = 6.31 nibal 



In 



/ 9/i„ 



V ^mbol 



-(35 - 250) MeV. 



2.06 
(3.4) 



To obtain the estimate for the ultrasoft constant 
terms, I interpreted as as(/^i)'*as(Mus), with 
(ii K, 2 GeV the scale that makes the logarithm / 



in (3.1) vanish, and varied /x„s from 0.7 GeV to 
2 GcV. The numerical estimate is highly sensitive 
to the PNRQCD cut-off scale ^us which must 
cancel in a complete order calculation. 
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The second estimate is based on the logarith- 
mically enhanced terms of order In as [^it |5^ : 



1730 



8l7r 

(75 - 100) MeV 



mbo;^ ln(l/as 



(3.5) 



Here I varied the scale in as\n{l/as) from 0.7 
to 2 GeV as above. This set of terms docs not 
depend on arbitrary cut-offs, but the logarithm 
is not large. Note that the coefRcient of the loga- 



rithm is not identical to the one in (3^), because 
there are logarithms unrelated to ultrasoft effects 
in the potentials of PNRQCD. 

Both corrections may not be small compared 



to the error estimate in ( 3.2 ) , but since they come 
with opposite sign and constitute only part of the 



al correction it is too early to revise (3^). How- 
ever, they illustrate that next-to-next-to-next-to- 
leading order effects may be large. 

3.2 Decay into a lepton pair 

The decay of a nS state into l^l^ measures the 
quarkonimn wave function ^'(0) at the origin. In 
turn, this parameter enters all those quarkonium 
decays which proceed through QQ annihilation. 
In terms of the quark electric charge eg, the fine 
structure constant a and the mass Mns of the 
quarkonium, the decay rate reads, for massless 
leptons: 



^^C(a.;/i)2|*(0)P(A.). 



(3.6) 



Here C(as;A*) is the short-distance coefficient of 
the non-relativistic vector current V'^cr^x, which 
is known to NNLO [|2[ ||, and ^-(0) is related 
to the NRQCD matrix element of the current. 
Note that the wave function at the origin is fac- 



torization scale dependent at NNLO. Eq. (3.6) 
neglects corrections from higher dimension oper- 
ators. They are incorporated in the numerical 
result below. 

The short-distance coefficient is poorly con- 
vergent in the MS factorization scheme implied 
by the threshold expansion. Numerically, at the 
NRQCD matching scale /i/i — mq, 



C{as]mQ) = 1 - 0.849 asimq) 
- (4.51 - QM2nf)as{mQf 



where Uf refers to all lighter flavours, approxi- 
mating them as massless. This provided the first 
hint that NNLO corrections to QQ systems are 
very large . 

The coefficient function is scheme-dependent 
and the large NNLO correction may be a scheme 
artefact. For perturbative QQ systems one can 
also compute |*(0)p(^) in PNRQCD perturba- 
tion theory (which implies treating the Coulomb 
interaction non-perturbatively). The NNLO 
correction to |^'(0)p(/i) has been obtained ana- 
lytically in and has been confirmed by [l6| . 
(Note, however, that the result is not presented 
explicitly in ||l^ . In Q a more complicated rep- 
resentation is given. According to ji^ the two 
representations agree numerically.) The decay 
width for T(IS') ^ l+l- to NNLO in PNRQCD 
perturbation theory, including the higher dimen- 
sion operators mentioned above, reads 



r(T(is') ^ i+r) 



32 
27 



2 2 / \; 



1 



1.99^-2.00) as{fi} + {2Ml^ + 3.261 



11.19 + 7.43 ln(m6/^)) as{fJ.f 



(3.8) 



where I is defined as in (3T). The logarithm 
of rrib/fi is partly related to the fact that the 
natural scale in the short-distance coefRcient is 
rrib and not /ii. However, a consistent treatment 
of all such logarithms requires renormalization 
group methods and has not yet been given. For 
/ii such that I — (with = 5 GeV), the factor 
in brackets is 



F(/ii) = 1 - 0.60+ 1.63. 



(3.9) 



As in (B.l) the series is not convergent at all. 
We have been able to eliminate such behaviour 
as unphysical for the quarkonium mass, but we 
cannot apply the same reasoning here. It there- 
fore seems that the leptonic decay width cannot 
be predicted reliably in perturbation theory. 



(3.7) 



Despite (3/7) bottomonium and charmonium 
decays may well be reliably predicted in NRQCD 
with the NRQCD matrix elements treated as non- 
perturbative parameters. To obtain a definite 
conclusion one would have to compute another 
quarkonium decay, such as rjb — > 77 or T(IS') 
light hadrons, to NNLO. The ratio of decay rates 
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is given by a factorization scheme independent 
ratio of short-distance coefficients. It may be 



that such ratios are more convergent than (3.7). 

A subset of next-to-next-to-next-to-leading 
order contributions to the T(15) leptonic decay 
width is known. The ultrasoft correction at order 
contributes 



-6.81a^ 



In 



0.777 



(-0.15) - (-^0.09). 



(3.10) 



to F in (U) 111. The numerical range is com- 
puted with the same prescription used for (3.4). 
The leading logarithmic contribution at order 
is a double logarithm. It arises as a product of 
an ultrasoft or potential logarithm and a loga- 
rithm related to the anomalous dimension of the 
current ip'^a'^x- The correction is [p2[: 



5F 



In 



212 



97r 

-(0.08- 0.30). 



ln(l/a,)2 



(3.11) 



The numerical range is computed with the same 
prescription used for ( [3.51 ). Again these correc- 
tions are not small, but their impact is less severe 
given the already large uncertainty of F{pi) in 



(3.9) 



In conclusion, higher order corrections al- 
most always turn out to be large. The T(15) 
mass may be useful to determine the b quark 
mass. A reliable prediction of absolute decay 
widths appears improbable, but it remains to be 
checked whether ratios behave better. This con- 
clusion may be frustrating. However, the numer- 
ical analysis of the quarkonium mass and wave 
function at the origin provides important insight 
into the structure of corrections to inclusive top 
and bottom quark pair production near thresh- 
old. 

4. Top quark pair production near 
threshold 

Top quark pair production is one of the major 
physics cases for a first linear e+e" collider and 
has been studied extensively in this context . 
The threshold behaviour of the cross section can 
be used to determine the top quark mass with 



great precision - provided the theoretical predic- 
tion is accordingly accurate. 

Toponium would be the perfect candidate 
for perturbative applications of non-relativistic 
QCD, but nature has provided another compli- 
cation. The electroweak decay width Ft w r(t 
bW) increases as mf. In the standard model 
Tt ~ 1.4 GeV, of the same order of magnitude 
as the ultrasoft scale of toponium jsj] . 

Suppose (for a moment) that the top quark 
is stable and neglect the axial- vector coupling to 
the Z boson, which is suppressed near threshold. 
Then the tt production cross section is obtained 
from the correlation function 

n^<>(g^) = ('Jm'?!^ - 9^5M'')n('7^) 
= t J d4a;e^«-(0|T(j,,(x)j,(0))|0), (4.1) 

where j'^{x) = [t'y'^t]{x) is the top quark vec- 
tor current and s = the centre-of-mass energy 
squared. Defining the usual i?-ratio R = cr^f/cro 
(ctq = 47rag„/(3s), where aem is the electromag- 
netic coupling at the scale 2'mt), the relation is 

Aire? 

R= '-{l + az)l^IV"is + ie), (4.2) 

s 

where et — 2/3 is the top quark electric charge 
in units of the positron charge and az accounts 
for the vector coupling to the Z boson. (For the 
remainder of this section, I set = for sim- 
plicity.) Only the spatial components of the cur- 
rents contribute up to NNLO. In the following 
mt refers to the top quark pole mass. 

At leading order in PNRQCD perturbation 
theory the heavy quark current two-point func- 
tion is given by the first diagram in figure ^. The 
tt pair is created by the current, interacts instan- 
taneously through the (leading order) Coulomb 
potential, and is annihilated by the current. Fig- 
ure ^ is actually misleading. Because tt produc- 
tion is short-distance compared to the toponium 
scale mtas, the pair is created and destroyed at 
the same space-time point! The leading order 
result is 



His) 



2rn? 



(4.3) 



where E = ^/s — 2mt. The Green function at 
the origin is ultraviolet divergent. In dimensional 
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II 



Figure 4: PNRQCD perturbative diagrams for tiie 
heavy quark current correlation function. At leading 
order the current generates a QQ pair which propa- 
gates with the Coulomb Green function (first line). 
Black bars denote insertions of potentials. The last 
diagram contains an ultrasoft gluon exchange. Both 
diagrams in the last line are beyond NNLO. 



regularization, with MS subtractions, one finds 

" 1 
, 2A 

1 

2 



1 , -AmtE 
^ In t: — 



(4.4) 



where A — 2as/ {3y/—E/mt) and 7^ is Euler's 
constant. The cross section requires only the 
discontinuity of 11 , which is scheme-independent. 
The Green function contains a continuum at E > 
and poles at £' < 0, from which the energy 
and wave function at the origin of toponium res- 
onances can be extracted. 

How does the top decay width affect this re- 
sult? For top quarks with energy E ~ mtv'^ ~ Tt 
and momentum p ^ mtv, we can approximate 
the quark propagator (in the potential region) 

1 



mt 



1 



[l + 0{v)]. (4.5) 



E + irt~p^/{2mt) 

As expected, the width is a leading order effect, 
but can be taken into account at this order by 
substituting E E = E + iVt, where Tt is the 
gauge-independent on-shell decay width. This 
gives the classic leading order result of Be- 
cause the Green function is evaluated off the real 
axis, the toponium poles are smeared out. For 
Ft « 1.4 GeV, we expect to see a broad remnant 
of the 15* pole, but all higher resonances overlap 
and form a continuum. 

Next-to-leading order corrections to this re- 
sult have been known for some time [p6f. Other 



properties of the production process, such as top 
quark momentum distributions and asymmetries 
generated by interference of vector and axial- vec- 
tor contributions, have been studied in some de- 
tail | |5^ , sometimes with non-perturbative mod- 
ifications of the heavy quark potential that have 
little justification in the theoretical framework 
described in earlier sections. The recent devel- 
opment concerns the calculation of NNLO cor- 
rections in a systematic non-relativistic approach 
m, |l|, H, 11, In PNRQCD Feynman di- 
agrams these corrections are given by the two 
diagrams in the second line of figure ^ with po- 
tentials up to NNLO. These diagrams correspond 
to integrals of the form 



d'^p (Pp' (Pqi dPq2 
(2^(2^(2^(2^ 



Gcip,qi;E) 



(4.6) 



and generalizations with more than one insertion 
of an interaction potential dV. Triple insertions 
of potentials and ultrasoft gluon exchange are 
higher order and neglected. The new calcula- 
tions either employ a combination of numerical 
and analytical methods ^ 38 or are fully 
analytical Q . The numerical solution of the 
Schrodinger equation contains higher order cor- 
rections, because an infinite number of insertions 
of the potentials that are kept in the equation is 
included. This could be an advantage if it could 
be argued that these higher order corrections are 
the dominant ones. 

The first NNLO calculations Q reported 
large corrections to the peak position and nor- 
malization of the remnant of the IS* toponium 
resonance. The correction to the peak position 
could be explained as an artefact of on-shell renor- 
malization of the top quark mass [^9| . Subse- 
quent calculations incorporated this observation 
|l|, H, m and used one or the other of the "al- 
ternative" mass definitions discussed in Sect. ^. 
The NNLO pair production cross section near 
threshold in the PS scheme is shown in the up- 
per panel of figure ^. For comparison the result 
in the pole mass scheme is shown in the lower 
panel - for the very last time! 

The generic features of figure ^ are easily un- 
derstood in terms of the results of Sect. |[ since 
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Figure 5: (a) [upper panel] : The normalized cross 
section (virtual photon contribution only) in LO 
(short-dashed), NLO (short-long-dashed) and NNLO 
(solid) as function of S = y/s — 2mt,ps(20 GeV) 
(PS scheme, yLf — 20GeV). Input parameters: 
mt,ps(20GeV) = = 175 GeV, Ft = 1.40 GeV, 
as{mz) ~ 0.118. The three curves for each case 
refer to fi — {15(upper); 30(central); 60(lower)} GeV. 
(b) [lower panel]: As in (a), but in the pole mass 
scheme. Hence E = -^s — 2mt. Other parameters as 
above with mi,Ps(20GeV) ^ m*. Plot from |ll 



the behaviour of successive perturbative approx- 
imations in the vicinity of the peak reflects es- 
sentiaUy the perturbative expansion of the topo- 
nium IS* mass and wave function at the origin. 
The difference in the shift of the peak position in 
the pole and PS scheme is a direct consequence 
of the improved convergence of the perturbative 
expansion of Mig. The stability of the peak po- 
sition in the PS scheme implies that the PS mass 
(but not the pole mass) can be determined accu- 
rately from the measurement of the cross section. 
The PS mass determined in this way can also be 
related more reliably to the top quark MS mass, 
which is probably the most useful reference pa- 
rameter. The numerics of table |l| adapted to the 
top quark gives, for given Wit = 165 GeV (and 



asim) = 0.1083): 

mt = [165.0 -I- 7.58 + 1.62 -I- 0.51 

0.24 (est.)] GeV (4.7) 
mt,ps(20 GeV) = [165.0 -I- 6.66 + 1.20 + 0.28 

0.08 (est.)] GeV, (4.8) 

where the numbers refer to successive terms in 
the perturbative expansion. The difference in 
convergence is significant on the scale of 0.1 GeV 
set by the projected statistical uncertainty on the 
mass measurement. 



However, the perturbative expansion (3 
for the wave function at the origin squared re- 
mains poorly convergent even for quarks as heavy 
as top. This leads to a large modification of the 
normalization of the cross section near the reso- 
nance peak at NNLO and also to a large renor- 
malization scale dependence since the NNLO cor- 
rection is proportional to a^. The recent calcu- 
lations |l^, m, agree qualitatively on the be- 
haviour of the peak position and normalization. 
When the Schrodinger equation is solved numer- 
ically, the scale dependence of the peak normal- 
ization appears to be smaller than in figure ^ 
| |3^ . It is an open question which of the two scale 
dependences provides a realistic estimate of the 
present theoretical uncertainty. 

A word of reservation applies to the term 
"NNLO" . All current NNLO calculations account 
for the width of the top quark by the replace- 
ment E E + iVt or a prescription of simi- 
lar parametric accuracy. Beyond a leading or- 
der treatment of (4.5), counting Tt rritv'^, the 
self-energy has to be matched to better accuracy. 
The correction terms relate to the off-shell self- 
energy and carry electroweak gauge-dependence. 
A complete NNLO result in the presence of a 
width that scales as above therefore includes elec- 
troweak vertex corrections as well as single res- 
onant backgrounds and non-factorizable correc- 
tions to the physical WWbb final state. Although 
some non-factorizable corrections are known near 



threshold |61 and away from threshold |62 , a 
systematic treatment of these complications has 
not been attempted yet. Strictly speaking, the 
concept of the tt cross section is not defined at 
NNLO and the problem has to be formulated in 
terms of a particular final state such as WWbb. 
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One may expect that single-resonant and non- 
factorizable corrections are 'structureless', that 
is, do not exhibit a pronounced resonance peak. 
In this case, they would add to the already exist- 
ing normalization uncertainty, but would affect 
little the top quark mass measurement. 

Besides the total cross section, top quark mo- 
mentum distributions js^, |3^ and vector-axial- 
vector interference [ |60| have been investigated. 
These quantities are even more delicate in the 
presence of a finite top quark width, which de- 
serves further investigation. Top quark produc- 
tion near threshold in 77 collisions was consid- 
ered in [Q. The theoretical accuracy is less in 
this case, because the two-loop short-distance co- 
efficient has not yet been calculated. 

How large could the corrections to figure || 
be? A particularly interesting set of corrections 
is again related to the ultrasoft scale. If we de- 
fine the ultrasoft scale as the scale where the log- 
arithm in (^^) vanishes, we obtain ~ 3 GeV 



for top quarks. The scale fix at which I in (3.1) 
vanishes is 32.6 GeV. I then repeat the estimates 



(3.4), (K§, (3.1C) and (3.11) for top quarks, us- 



mg agam the resuhs of |2lJ |l], ||, and varying 
Hus between 2 GeV and 5 GeV. I obtain 



-(140- 300)MeV 
0.01-0.05 



for the ultrasoft correction and 



spin 



(150- 160) MeV 
-0.07 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



for the leading logarithmic NNNLO correction. 
These numbers are relevant to tt production in 
the vicinity of the resonance peak. The correc- 
tions to the peak position are again not small 
compared to the residual uncertainty of about 
200 MeV in figure |^, but the two corrections are 
of opposite sign and a conclusive estimate re- 
quires further NNNLO terms. In addition, one 
may worry about finite width effects in ultrasoft 
contributions. 

To conclude this section, let me make the 
following remark. It is often said that the width 
of the top quark screens non-perturbative QCD 
effects and makes the threshold cross section cal- 
culable in perturbative QCD. It is true that the 



width smears the toponium resonances and con- 
verts the observed cross section into a smooth 
excitation curve. However, it is not true that 
the top quark width screens non-perturbative ef- 
fects any more than the existence of a pertur- 
bative ultrasoft scale mtol ~ Ft already does. 
Even for stable top quarks the production cross 
section near threshold (averaged over an inter- 
val several times Aqcd) is perturbatively calcu- 
lable, as would be the toponium resonances and 
their decays. Interestingly, perturbative resum- 
mation with power counting v ^ Ug seems to 
make sense even as w — s- at fixed as ■ In particu- 
lar, the cross section directly at threshold seems 
to be infrared safe in perturbation theory. As 
u ^ the arguments of the coupling constants 
as(mQv) and as(mQW^) freeze at values of or- 
der niQas and raqol, respectively, and do not 
tend to zero. This has been checked explicitly 
by investigating the logarithms up to order 
(NNLO). 

5. Sum rules, the h quark mass 

As a final application we return to the b quark 
mass. There are legitimate doubts concerning 
the reliability of ( |3.3D given that T(IS') is hardly 
a truly perturbative onium. One can by-pass this 
problem by considering averages over the bottom 
pair production cross section rather than exclu- 
sive resonances. This leads us to consider sum 
rules lei 



M„/(10GeV) 



2n 



127r2 d" 



n! rf(g2)' 



■n(g^)i 



(5.1) 



which equate an experimental average of the cross 
section to the perturbatively computable deriva- 
tives ("moments") of the bottom vector current 
correlation function. 

The parameters of the lowest T{nS) reso- 
nances are well measured and one chooses n large 
enough that the experimental uncertainty on the 
continuum cross section is small compared to the 
theoretical uncertainty. This occurs for n > 6. 
For such moments the ordinary perturbative ex- 
pansion of the moments breaks down and has 
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to be replaced by non-relativistic resummation 
and PNRQCD perturbation theory. Leading or- 
der and next-to-leading order analyses of the sum 
rule with non-relativistic resummation have been 
performed in ||6^ and Q , respectively. After in- 
tegration over s the expansion of Rbtis) in Us/v 
turns into an expansion in ag^/n. Non-relativistic 
resummation sums these terms to all orders. For 
moments the scale uibV turns into 2mf, / -^/n; the 
ultrasoft scale mf,w^ is rrib/n. The requirement 
that the ultrasoft scale is perturbative translates 
into n < 10. Larger moments are often used in 
the literature. This introduces a systematic un- 
certainty which is difhcult to quantify at NNLO, 
since at this order ultrasoft corrections are not 
included. 

Several NNLO calculations have been com- 
pleted recently |||, ||, |5^, |4|, the calculation 
being almost identical to that of top pair pro- 
duction near threshold. The later publications 
||50| , H^, abandoned the idea of determining 
the pole mass and usually extract the bottom 
MS mass mi,. This is done by extracting the PS 
or IS or kinetic mass from the sum rule, which 
is then converted into rub. Although the differ- 
ent groups compute the same quantity, there are 
differences in the implementation of the resum- 
mation which are formally beyond NNLO. These 
concern (a) whether the short-distance coefficient 



(3.7) is kept as an overall factor or multiplied out 



to NNLO; (b) whether the integral over s in (5.1) 
is done exactly or in a non-relativistic approxima- 
tion; (c) whether the energy denominator of the 
full Green function is expanded up to NNLO or 
whether the exact NNLO energy levels are kept. 
These differences in implementation can shift the 
value of TOfc extracted from the sum rule by up 
to 100 MeV. 

Further differences arise in the choice of mo- 
ments, renormalization scale, at which the sum 
rule is evaluated, and the analysis strategy. 

Refs. 1^ perform essentially a single- 
moment analysis, based on the fact that the the- 
oretical error is highly correlated between differ- 
ent moments. It is then checked that varying the 
moment gives a negligible difference. One then 
finds a significant renormalization scale uncer- 
tainty, which can be traced back to the badly be- 




4 6 8 
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Figure 6: (a) [upper panel]: The value of 
mi,^Ps(2GeV) obtained from the 10th moment as a 
function of the renormalization scale in NLO and 
NNLO and for as{mz) = 0.118. The dark region 
specifies the variation due to the experimental error 
on the moment. The middle line marks the scale 
fi„ = 2mi,/y'n, the two outer lines determine the 
scale variation from which the theoretical error is 
computed. Plot from jisj. (b) [lower panel]: The 
value of mb,ps(2 GeV) obtained from a linear combi- 
nation of moments as a function of the renormaliza- 
tion scale in NNLO. M„ is defined with the normal- 



ization of (5.1) 



haved expansion (3^). A typical result is shown 



in the upper graph of figure |g. The results 

mf,,ps(2 GeV) = (4.60 ±0.11) GeV @ (5.2) 
m6,kin(l GeV) = (4.56 ± 0.06) GeV [|o) (5.3) 

differ by about 80 MeV, when related to each 
other according to table |l|, but are consistent 
with each other within implementational differ- 
ences. The larger error on the first result follows 
from a larger renormalization scale variation. 

The analysis of |Q is different, because it 
uses linear combinations of moments. The linear 
combination chosen in [Q is less sensitive to a 
variation of the renormalization scale, while re- 
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Rcfs. 


mb{mb) 


nib 


Remarks 


@ 

0, |6l 

ii 


4.20 ±0.10 
4.19 ± 0.06 
4.26 ±0.10 

4.21 ± 0.11 


4.97±0.17 
4.80 ±0.06 


Sum rules (via mb_kin(l GeV)) 

Sum rules (via mb is) 
Sum rules (via TOf,,ps(2 GeV)) 
Sum rules (via nib, 2-loop) 


ii 
il 
11 


4.44 ±0.04 
4.24 ±0.09 
4.21 ±0.07 


5.04 ±0.09 


T(lS')mass (via nib, 2-loop) 
T(IS') mass (via TOf,^ps(2 GeV)) 
T (IS*) mass (direct) 



Table 3: Bottom quark MS and pole mass values (in GeV) obtained from NNLO sum rule or T(15') mass 
calculations. 



taining sensitivity to nib. This is illustrated in 
the lower panel of figure ^ in the PS scheme. In 
this way, Q obtains 



mb,is = (4.71 ±0.03) GeV, 



(5.4) 



which is close to ( |5.3| ) after conversion. This re- 
sult depends crucially on combining moments at 
identical renormalization scales rather than their 
"natural" scale 2nib/\/n, and on discarding pos- 
sible multiple solutions (as the upper one in fig- 
ure ^). In my opinion the error quoted in (5.4) 
should be understood as an error that follows un- 
der the specific assumptions of the analysis strat- 
egy. It can hardly be considered as a realistic 
estimate of the total theoretical error, given the 
differences that can arise in different implemen- 
tations of the theoretical moments and given the 



size of ultrasoft effects discussed in Sect. 3.1 



The results quoted above can be converted 
into the bottom MS mass using table |l| A sum- 
mary of NNLO results from sum rules and, for 
comparison, the T(15') mass, is compiled in ta- 
ble ^, where I quote the number given by the 
authors. This number may differ from the one 
of table [l], because not always is a four-loop re- 
lated to TOfc used as appropriate to a NNLO sum 
rule calculation (cf. the discussion after (|3.3|)). 
This difference is small when mass definitions 
with convergent relations to fub are used, see ta- 
ble ^, but affects Wib when computed from the 
pole mass via a 2-loop relation. For this rea- 
son the results for nib from |49| in table || 
should in fact be decreased by about 200 MeV. 
This makes consistent with the other rrib de- 
terminations, but puts the result of [}49| off by 



200 MeV. Comparison of pole mass results shows 
that this discrepancy is already present in the 
pole mass value before conversion to mb- The 
small pole mass value of |^ is a consequence of 
the fact that the sum rule is evaluated at high 
renormalization scales. Fig. |^ shows that this 
leads to reduction of nib- 

My "best" estimate for the b quark MS mass 
is a (potentially biased) combination of the re- 
sults of [|o[ H, Hi. It is remarkable that the 
sum rule result is consistent with the T(IS') re- 
sults despite the fact that the systematics of non- 
perturbative effects and scale dependence is dif- 
ferent. Nevertheless, the results of table ^ are 
not independent and may be affected by com- 
mon, unidentified theoretical uncertainties. This 
being said, my preferred value is 

mbimb) = (4.23 ± 0.08) GeV. (5.5) 

This is in beautiful agreement with nibijiib) = 
4.26±0.07GeV Q obtained from lattice HQET. 
This uses the B meson mass, a lattice calcu- 
lation of the (properly defined) binding energy 
of the B meson in the unquenched, two- flavour 
approximation to heavy quark effective theory, 
and a two-loop perturbative matching to the MS 
scheme. To our knowledge, this is the only other 
NNLO determination of the MS mass besides the 
sum rule calculations mentioned above (which, 
in fact, are N^LO as far as rfib is concerned). In 
my opinion using an error smaller than 80 MeV 
on fubirnb) or any of the "alternative" masses 
cannot be justified at present. Using smaller er- 
rors in B physics observables may have danger- 
ous consequences for consistency checks of the 
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CKM model of CP violation. 

6. Concluding remarks 

The theory of pcrturbative onium systems has 
imdergnne a,n eyeiting tra.nsitinn from potential 



jdels with arbitrar}^ cut-offs and poorl}^ imder- 



stood accuracy to a systematic effective theory 
description. This transition can be compared in 
significance with the development of heavy quark 
effective theory for B and D mesons. Unfortu- 
nately, nature has not been kind to us, leaving us 
with systems which are barely pcrturbative (bot- 



tomonium) or extremely short-lived (toponium). 



Along with this development pcrturbative ex- 
pansion tools have been invented and the most 
basic quantities are now computed to ncxt-to- 
next-to-leading order. These calculations sharp- 
ened our understanding of heavy quark mass re- 
normalization, but large corrections remain in 
most cases. They confront us with the challenge 
of a complete next-to-next-to-next-to-leading or- 
der calculation. With so many tools at hand, 
progress is surely expected. 

Acknowledgments 

I would like to thank the organizers of HF8 for 
making this a focused and enjoyable meeting. 
Thanks to A. Signer and V.A. Smirnov for collab- 
oration on subjects related to this review and to 
A.H. Hoang and A. Pineda for discussions. This 
work was supported in part by the EU Fourth 
Framework Programme 'Training and Mobility 
of Researchers', Network 'Quantum Chromody- 
namics and the Deep Structure of Elementary 
Particles', contract FMRX-CT98-0194 (DG 12 - 
MIHT). 

References 

[1] E. Eichten and B. Hill, Phys. Lett. B234 
(1990) 511; B. Grinstein, Nud. Pfiys. B339 
(1990) 253; H. Georgi, Phys. Lett. B240 (1990) 
447. 

[2] N. Isgur and M.B. Wise, Phys. Lett. B232 
(1989) 113; B237 (1990) 527. 

[3] W.E. Caswell and G.P. Lepage, Phys. Lett. 
B167 (1986) 437. 



[4] B.A. Thacker and G.P. Lepage, Phys. Rev. 
D43 (1991) 196; G.P. Lepage et al, Phys. Rev. 
D46 (1992) 4052 

[5] G.T. Bodwin, E. Braaten and G.P. Lepage, 
Phys. Rev. D51 (1995) 1125; D55 (1997) 5853 
(E); LZ. Rothstein, these proceedings hep- 



[9: 

[10 

[11 
[12: 

[13 
[14 

[15: 
[16: 

[17 

[18 
[19 

[20 

[21 
[22: 



ph/991127e| 



A.V. Manohar, Phys. Rev. D56 (1997) 230. 

M. Beneke, Non-relativistic effective theory for 
quarkonium production in hadron collisions in: 
The strong interaction, from hadrons to par- 
tons, p. 549, Stanford 1996 ||hep-ph/9703429 |. 



P. Labelle, Phys. Rev. D58 (1998) 093013 hep- 



ph/9608491| 



M. Luke and A.V. Manohar, Phys. Rev. D55 

(1997) 4129. 

B. Grinstein and LZ. Rothstein, Phys. Rev. 
D57 (1998) 78. 

M. Luke and M.J. Savage, Phys. Rev. D57 

(1998) 413. 

M. Beneke and V.A. Smirnov, Nucl. Phys. 
B522 (1998) 321. 

H.W. Griefihammer, Phys. Rev. D58 (1998) 
094027; Power counting and beta function in 
NRQCD ]hep-ph/9810"23^ . 

A. Pineda and J. Soto, Nucl. Phys. Proc. Suppl. 
64 (1998) 428 |hep-ph/970748l[ ; Phys. Rev. 
D59 (1999) 016005. 

M. Beneke, New results on heavy quarks near 



threshold [ |hep-ph/9806429 |. 

M. Beneke, A. Signer and V.A. Smirnov, Phys. 
Lett. B454 (1999) 137. 

M.E. Luke, A.V. Manohar and LZ. Roth- 
stein, Renormalization group scaling in non- 



relativistic QCD |hep-ph/9910209 |. 

Y. Schroder, Phys. Lett. B447 (1999) 321. 

T. Appelquist, M. Dine and I.J. Muzinich, 
Phys. Lett. 69B (1977) 231; Phys. Rev. D17 
(1978) 2074. 

N. Brambilla, A. Pineda, J. Soto and A. Vairo, 
Phys. Rev. D60 (1999) 091502; Potential 
NRQCD: An effective theory for heavy quarko- 
nium [^ep|ph7990724^. 

N. Brambilla, A. Pineda, J. Soto and A. Vairo, 
The heavy quarkonium spectrum at order 



mallna^ [hep-ph/991023^ 



M. Beneke, A. Signer and V.A. Smirnov, Phys. 
Rev. Lett. 80 (1998) 2535. 



20 



Heavy Flavours 8, Southampton, UK, 1999 



M. Beneke 



A. Czarnecki and K. Melnikov, Phys. Rev. Lett. 
80 (1998) 2531. 

H. Leutwyler, Ptiys. Lett. B98 (1981) 447; 
M.B. Voloshin, Sov. J. Nucl. Phys. 36(1) 
(1982) 143; A. Pineda, Nucl. Phys. B494 
(1997) 213. 

A.S. Kronfeld, Phys. Rev. D58 (1998) 051501. 

I. I. Bigi and N.G. Uraltsev, Phys. Lett. B321 
(1994) 412. 

M. Beneke and V.M. Braun, Nucl. Phys. B426 

(1994) 301; M. Beneke, Phys. Lett. B344 

(1995) 341. 

I.I. Bigi, M.A. Shifman, N.G. Uraltsev and 
A.I. Vainshtein, Phys. Rev. D50 (1994) 2234. 

M. Beneke, Phys. Lett. B434 (1998) 115. 

A.H. Hoang, M.C. Smith, T. Stelzer and 
S. Willenbrock, Phys. Rev. D59 (1999) 
114014; A. Pineda, Heavy quarkomum and 
nonrelativistic effective theories, PhD thesis, 
Barcelona (1998). 

M.C. Smith and S.S. Willenbrock, Phys. Rev. 
Lett. 79 (1997) 3825. 

N. Gray, D.J. Broadhurst, W. Grafe and 
K. Schilcher, Z. Phys. C48 (1990) 673. 

K.G. Chetyrkin and M. Steinhauser, Short dis- 
tance mass of a heavy quark at order [hep 



ph/9907509| 



M. Beneke and V.M. Braun, Phys. Lett. 
B348 (1995) 513; P. Ball, M. Beneke and 
V.M. Braun, Nucl. Phys. B452 (1995) 563. 

See, for example, the review: M. Beneke, Phys. 
Rep. 317 (1999) 1. 

U. Aglietti and Z. Ligeti, Phys. Lett. B364 
(1995) 75. 

A.H. Hoang, Z. Ligeti and A.V. Manohar, 
Phys. Rev. Lett. 82 (1999) 277; Phys. Rev. 
D59 (1999) 074017. 

A.H. Hoang and T. Teubner, Phys. Rev. D60 
(1999) 114027. 

A. Pineda and F.J. Yndurain, Phys. Rev. D58 
(1998) 094022. 

I.I. Bigi, M.A. Shifman and N.G. Uraltsev, 
Ann. Rev. Nucl. Part. Set. 47 (1997) 591. 

A. Czarnecki, K. Melnikov and N. Uraltsev, 
Phys. Rev. Lett. 80 (1998) 3189. 

A.H. Hoang, IS and MS bottom quark m asses 
from Upsilon sum rules \ liep-ph/990555C | . 



[43; 

[44; 
[45; 

[46; 
[47 
[48; 
[49 

[5o; 

[51 
[52; 

[53; 

[54; 
[55; 

[56; 



M. Beneke and A. Signer, The bottom MS 
quark mass from sum rules at next-to-next-to- 
order [|hep-ph/990647e |. 



N.G. Uraltsev, Int. J. Mod. Phys. All (1996) 
515. 

M. Beneke, V.M. Braun and V.I. Zakharov, 
Phys. Rev. Lett. 73 (1994) 3058; A. Sinkovics, 
R. Akhoury and V.I. Zakharov, Phys. Rev. 
D58 (1998) 114025. 

P. Ball, M. Beneke and V.M. Braun, Phys. Rev. 
D52 (1995) 3929. 

T. van Ritbergen, Phys. Lett. B454 (1999) 
353. 

I.I. Bigi, M.A. Shifman, N.G. Uraltsev and 
A.I. Vainshtein, Phys. Rev. D56 (1997) 4017. 

A. A. Renin and A. A. Pivovarov, Phys. Lett. 
B435 (1998) 413; Nucl. Phys. B549 (1999) 
217. 

K. Melnikov and A. Yelkhovsky, Phys. Rev. 
D59 (1999) 114009. 

B. A. Kniehl and A. Renin, Ultrasoft effects in 
heavy-quarkomum physics | hep-ph/9907489|. 

B.A. Kniehl and A. Renin, Order aJln^(l/Qs) 
corrections to heavy- quarkonium creation and 



hilation hep-ph/9911414 



ECFA/DESY LC Physics Working Group 
(E. Accomando et al), Phys. Rep. 299 (1998) 
1. 

I.I. Bigi et al, Phys. Lett. B181 (1986) 157. 

V.S. Fadin and V.A. Khoze, Pts'ma Zh. Eksp. 
Teor. Fiz. 46, 417 (1987) [JETP Lett. 46, 525 
(1987)]; Yad. Fiz. 48, 487 (1988) [Sov. J. Nucl. 
Phys. 48(2), 309 (1988)]. 

M.J. Strassler and M.E. Peskin, Phys. Rev. 
D43 (1991) 1500. 



[57] W. Kwong, Phys. Rev. D43 (1991) 1488; 
M. Jezabek, J.H. Kiihn and T. Teubner, Z. 
Phys. C56 (1992) 653; H. Murayama and 
Y. Sumino, Phys. Rev. D47 (1993) 82; M. Jez- 
abek and T. Teubner, Z. Phys. C59 (1993) 
669; Y. Sumino, K. Fujii, K. Hagiwara, H. Mu- 
rayama and C.-K. Ng, Phys. Rev. D47 (1993) 
56; K. Fujii, T. Matsui and Y. Sumino, Phys. 
Rev. D50 (1994) 4341; R. Harlander, M. Jez- 
abek, J.H. Kiihn and T. Teubner, Phys. Lett. 
B346 (1995) 137; M. Jezabek et al, Phys. Rev. 
D58 (1998) 014006. 



21 



Heavy Flavours 8, Southampton, UK, 1999 



M. Beneke 



[58] A.H. Hoang and T. Teubner, Phys. Rev. D58 
(1998) 114023; K. Melnikov and A. Yelkhovsky, 
Nucl. Phys. B528 (1998) 59; O. Yakovlev, 
Phys. Lett. B457 (1999) 170. 

[59] T. Nagano, A. Ota and Y. Sumino, Phys. Rev. 
60 (1999) 114014. 

[60] A. A. Penin and A. A. Pivovarov, Analytical re- 
sults for e^e~ — > ti and 77 tt observables 
near the threshold up to the next-to-next-to- 
leading order of NRQCD [ |hep-ph/990427S ]. 

[61] M. Peter and Y. Sumino, Phys. Rev. D57 
(1998) 6912. 

[62] K. Melnikov and O. Yakovlev, Phys. Lett. 
B324 (1994) 217; V.S. Fadin, V.A. Khoze and 
A.D. Martin, Phys. Rev. D49 (1994) 2247; 
W. Beenakker, F.A. Berends, A. P. Chapovsky, 
Phys. Lett. B454 (1999) 129. 

[63] A. A. Penin and A. A. Pivovarov, Nucl. Phys. 
B550 (1999) 375. 

[64] V.A. Novikov et al, Phys. Rev. Lett. 38 (1977) 
626; 38 (1977) 791 (E); Phys. Rep. 41 (1978) 
1. 

[65] M.B. Voloshin and Yu.M. Zaitsev, Usp. Fiz. 
Nauk. 152 (1987) 361 [Sov. Phys. Usp. 30(7) 
(1987) 553]. 

[66] M.B. Voloshin, Int. J. Mod. Phys. AlO (1995) 
2865. 

[67] A.H. Hoang, Phys. Rev. D59 (1999) 014039. 
[68] A.H. Hoang, Heavy quark masses from the 



QQ threshold and the upsilon expansion [hep- 
pli/9909356t. 



[69] M. Beneke, A. Signer and V.A. Smirnov, A 
two-loop application of the threshold expansion: 
The bottom quark mass from b anti-b produc- 
tion [ ^ep-pli/9906476 1. 

[70] V. Gimenez, L. Giusti, F. Rapuano and 
G. Martinelli, NNLO unquenched calculation 
of the b quark mass [iep-lat/9909138]. 



22 



